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Using tlie density matrix renormalization group teclmique, 
we calculate numerically the low energy excitation spectrum 
and magnetization curve of the spin-1 antiferromagnetic chain 
in a staggered magnetic field, which is expected to describe 
the physics of R2BaNiOe,{R ^ Y) family below the Neel tem- 
perature of the magnetic rare-earth [R) sublattice. These re- 
sults are valid in the entire range of the staggered field, and 
agree with those given by the non-linear a model study for 
small fields, but differ from the latter for large fields. They are 
consistent with the available experimental data. The correla- 
tion functions for this model are also calculated. The trans- 
verse correlations display the anticipated exponential decay 
with shorter correlation length, while the longitudinal corre- 
lations show explicitly the induced staggered magnetization. 

PACS Numbers: 75.10.Jm, 75.40. Mg 



The quasi-one-dimensional magnets have been the fo- 
cus of analytic, numerical and experimental studies since 
Haldane pointed out the difference between the integer 
spin Heisenberg antiferromaeiietic (AF) chains and the 
half-integer chains in 1983. M By mapping the Heisea- 
berg spin chains onto the 0(3) nonlinear cr-model, 
he conjectured that the low-energy excitation spectrum 
displays a finite gap for the integer spin systems while 
it is gapless for half-integer spin chains. This conjec- 
ture has been verified by later experiments on quasi- 
one-dimensional spin-1 materials such as NENP and 
Y2Baffi05 which show clear evidence of the Haldane 
gap. B Nowadays, the pure one-dimensional Haldane 
systems are fairly well understood, and a reliable esti- 
mate for the Haldane gap A = 0.41048(2)J for spin-1 
chains has been obtained by both dcnsi|t^ matrix renor- 
malization group (DMRG) calculation B and finite size 
exact diagonalization. M 

More recent developments on the Haldane systems con- 
cern various effects of external perturbations: doping 
with magnetic or non-magnetic impurities and apply- 
ing external magnetic field. The impurity doping may 
introduce bound states within the Haldane gap, ID while 
applying uniform external magnetic field splits the degen- 
erate HaldaBC triplet state into transverse and longitudi- 
nal modes. B The longitudinal mode becomes softened 
upon increase of the magnetic field, and at a critical field 



He the system enters a new phase with long-range AF 
order. Of course, a staggered applied magnetic field is 
even more interesting which would induce non- vanishing 
staggered magnetization and affect the Haldane gap ex- 
citation spectrum, but such a staggered field cannot be 
materialized by an external source. 

Most recently, a series of experiments performed on the 
family of quasi-one-dimensional materials with a general 
formula R2BaNi05, Wt3 where R is one of the mag- 
netic rare-earth elements substituting fully or partially 
Y (for brevity we denote this replacement by i? ^ Y), 
have made it possible to study the effect of the stag- 
gered magnetic field on the Haldane systems in detail. 
All members of this family contain spin-1 Ni"^^ linear 
chains and the in-chain AF exchange coupling is rather 
strong. (The detailed structure of this family of com- 
pounds is described in Ref. ||ll|). The reference com- 
pound Y2BaNiO^ is found to be highly ane-dimensional 
with negligible interchain interactions, o and no mag- 
netic order has_been observed so far even at very low 
temperatures. Hence it is believed to be an almost 
ideal example of the Haldane-gap system. Other mem- 
bers have magnetic R^~^ ions in addition to the spin-1 
Ni^^ ions. These ions are positioned between two neigh- 
boring Ni chains, weakly coupled to the iVi^+ ions and 
the coupling between themselves is also very weak. Nev- 
ertheless, these magnetic ions are AF ordered below cer- 
tain Neel temperature Tn- These ions do not affect the 
Ni chains substantially above Tjsi , keeping their Haldane 
features untouched, but the 3D AF ordered R^^ sublat- 
tice below Tiq has dramatical effects on these chains, im- 
posing a staggered magnetic field. The neutron scatter- 
ing experiments on powder samples and small. siac-.single- 



crystals of Nd2BaNiOz and Pr2BaNiO^ 



show 



an increase of the energy gap below the Neel tempera- 
ture. 

We assume that these chains can be still considered 
one-dimensional, being put in a staggered field created by 
3D ordered ions at low temperatures. The Hamilto- 
nian can be then writen as : 



(1) 



where J is the exchange constant (to be taken as en- 
ergy unit, i.e., J =1 ). The dimensionless staggered field 



1 



h = gS^BH-ir/ J with iJ^r as the physical staggered field, 
which, in turn, is proportional to the R sublattice mag- 
netization M/j 
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H„ = uMr. 



(2) 



g=2 is the theoretically predicted gyromagnetic ratio of 
the Ni ion. The direction of the staggered field has been 
chosen as the z axis. This Hamiltopian. has been con- 
sidered using the mean-field theory O'Ej' as well asJajy 
mapping onto the 0(3) nonlinear a model (NLSM). [E2l 

In this Communication we use the DMRG technique to 
calculate the low energy excitation spectrum and magne- 
tization curve of the Hamiltonian Eq. ([|) . The obtained 
field dependence of the gap and the staggered magneti- 
zation is consistent with the experimental results. These 
results are valid in the entire range of the staggered field 
and recover those given by non-linear a model for small 
fields but differ from the latter for large fields. More- 
over, we calculate the spin-spin correlation functions for 
this model. The transverse correlations display an ex- 
ponential decay as anticipated for the spin-1 AF chain, 
with a shorter correlation length, while the longitudinal 
correlations show explicitly the induced staggered mag- 
netization, nnn 

We follow the standard DMRG algorithm mM5 to 
calculate the low-energy excitations of the Hamiltonian 
(|l|), adopting the periodic boundary conditions (PBC). 
We use the infinite-chain algorithm up to chain length 
N—&Q and keep as many as 400 optimized states during 
each sweep. The largest truncation errors are of the order 
of 10~^ for smaller h, while for bigger h, these errors are 
as small as 10"^'^, which means our results are even more 
reliable for bigger h. 

The numerical results for the change of the lowest ex- 
citation energies (the Haldane gap) of Hamiltonain (|^) 
A — Ao are presented in Fig. |l| as functions of the dimen- 
sionless staggered magnetic field h. In the absence of this 
field the longitudinal (A^) and transverse (At) modes 
are degenerate, forming the Haldane triplet. For non- 
zero h, these modes will split with respect to each other. 
Both of them will increase with the staggered magnetic 
field, while the longitudinal gap increases faster than the 
transverse one. For small staggered fields, the increase 
of the longitudinal gap will be nearly three times faster 
then the transverse ones, while for larger staggered fields, 
this ratio will decrease, and is is approximately two for 
the largest staggered field we considered. 
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FIG. 1. The DMRG results for the transverse (solid circle) 
and longitudinal (empty circle) energy gaps as functions of the 
staggered magnetic field for spin-1 chain. 

The staggered magnetic moment of the system with 
chain length N is defined as 



M^{N) 



-Yi- 



(3) 



the longest chain in our calculation being N=60. 
Then the staggered magnetic moment for infinite long 
chain can be obtained by extrapolating Af^r = 
limjv^oo Mt^IN). Obviously, this quantity is a function of 
the staggered field, and our numerical results are shown 
in Fig. |2[ Considering Eq. (||), this figure is nothing 
but the relation between the magnetization of the Ni 
sublattice and the R sublattice, and it is qualitatively 
in agreement with the exprimental data in Fig. ^ of 
Ref. jl6| . Our numerical results show that in small stag- 
gered fields, the magnetization change linearly with the 
increase of the field, so we can easily extract the "zero 
field" staggered magnetic susceptibility x''*''(0)=18.50/J. 
This value fully agrees with the results obtained from 
the transfer-matrix renormalizaion group c3 , Quantum 
Monte Carlo B and the NLSM calculations. H We fit 
our results using the following function: 

= a arctan(fe /i) + (1 - ^a) tanh(c /i'^) (4) 

with a=0.412, 6==38.106, c==1.195, d=0.621. The fittig 
line is alsa,shown in Fig. To compare with the NLSM 
results 62! in detail, in Fig. ||, we use their analytic 
relation 
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{0)h = M^(l + 1.56A/2 + 2.4Af^ -I- 3.27M^) 



with their value x^'^HO) — 18. 7/ J as a reference. We see 
clearly that the analytic expression is very good for small 
staggered fields, while for larger staggered fields it devi- 
ates from our numerical results significantly. We have 
also calculated the magnetic moment for large h which 
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is not shown in Fig. Our result indicates that the 
moment should saturate in large enough staggered mag- 
netic field for both isotropic and (single-ion) anisotropic 
cases. The unsaturated moment at zero-temperature ob- 
served so far in various experiments tells us that the in- 
duced staggered magnetic field on Ni chains is, proba- 
bly, not large enough yet. This issue was also discussed 
earlier. H 
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FIG. 2. The staggered magnetization curve for spin-1 
chain. The staggered magnetic field is proportional to the 
magnetization of the R sublattice. The numerical results 
(solid circle) are fitted by a function with four parameters 
(solid line) (see Eq. (||)), with a=0.4f2, 6=38.106, c=1.195, 
d= 0.621, respectively; the analytic results given by non-linear 
a model are also presented (dotted line). 

From the above results, we obtain the values of the 
transverse as well as the longitudinal gap as functions of 
the magnetic moment of the Ni sites, which can be com- 
pared directly with the analytic NLSM result (see Fig. 
^). These xeaults are also consistent with the experimen- 
tal data. As for the comparison with the NLSM 
treatment we see again that in small staggered fields, 
the analytic and numerical results are in good agree- 
ment with each other, while for larger staggered fields, 
the analytic results deviate significantly from the numer- 
ical simulations. Both longitudinal and transverse gaps 
increase faster in simulations than the NLSM predicts. 
Since for our DMRG calculations, larger is the staggered 
field, more reliable are the results, so the disagreement of 
magnetization moment and the gaps between these two 
approaches raises a question whether the NLSM mapping 
is valid or not for large staggered magnetic fields. 
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FIG. 3. The numerical results for the transverse gaps 
(solid circle) and the longitudinal gaps (empty circle) ver- 
sus the magnetic moment on Ni sites; the NLSM results are 
shown ( solid line for transverse gaps and dotted line for lon- 
gitudinal gaps) in the same figure. 

Besides calculating the low-energy spectrum of the 
one-dimensional systems, the DMRG method also pro- 
vides a direct and simple way to calculate the spin-spin 
correlation functions which can shed some further light 
on the nature of the system under study. For a pure Hal- 
dane system, the correlation dacays exponentially, fol- 
lowing 



(^0 



x{z) gx{z)^ 



(5) 



where ^=6.03 is t^e correlation length obtained by the 



numerical study, and A is a constant. When a stag- 
gered field is applied, the AF long-range order will be in- 
duced along the z-direction, so {S^Sf) will not decay any 
more. However, {S§S[) — {Sq){Si) will still decay expo- 
nentially. Of course, the transverse correlations (SgSf) 
decay exponentially as before, but with modified expo- 
nents. In Fig. ||(a) and Fig. ^(b), the functions Cxx{l) = 
HVlKSSSm and a.(0 = HVl\{S^S[) - {S§){S[)\) 
are shown for different staggered fields. We see that both 
functions decay exponentially following Eq. (||) and the 
correlation lengths (,xx and £,zz decrease with increasing 
staggered field. As h increases, the reduced longitudi- 
nal correlations decay much faster. We have extracted 
the correlation lengths and extrapolated them to infinite 
chain length by considering the chain length dependence 
of ^xx and ^zz ■ We find that At and ^ Al , 

and both resuhs can be fitted by 0.402 * A (Fig. 

|), which coincides exactly with ^ ~ 6.03 ip)3tained for 
the isotropic spin chain {h = 0, Aq — 0.41). Q This is an 
independent check of the self-consistency in our calcula- 
tions. 
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In conclusion, by considering a model hamiltonian 
which describes the physics of a family of mixed-spin ma- 
terials in the temperature range below Tm of the mag- 
netic rare-earth sublattice, we have calculated numeri- 
cally the energy gap and the staggered magnetic moment 
as functions of the staggered magnetic field created by the 
AF long range order. The obtained results are consistent 
with the experimental data qualitatively. Our numerical 
results are also compared with the analytic considera- 
tions based on the non-linear a model. The comparison 
shows that the NLSM results are good for small stag- 
gered fields, while they deviate from the numerical simu- 
lations for larger staggered fields. After submitting this 
paper we saw a new report on polarized neutron studyj-of 
longitudinal Haldane-gap excitations in Nd2BaNi05. E3 
which show somewhat different behavior than expected 
from the theory. The reason of this discrepancy has to 
be understood. 
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FIG. 4. Spin-spin correlation functions for different val- 
ues of staggered magnetic field, where I is the distance 
between the two spins. Both correlations (SgSf) and 
(SoSi) — {So){Si} decay exponentially, but the latter decays 
faster. 
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FIG. 5. The inverse correlation lengths and ^^^^ for 
infinite chain length vs. the transverse and longitudinal gaps, 
respectively. The solid line is the fiting (^^^ — 0.402 * A. 
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